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Abstract 



We develop a linear theory of electron transport for a system of two iden- 
tical quantum wires in a wide range of the wire length L, unifying both the 
ballistic and diffusive transport regimes. The microscopic model, involving the 
interaction of electrons with each other and with bulk acoustical phonons al- 
lows a reduction of the quantum kinetic equation to a set of coupled equations 
for the local chemical potentials for forward- and backward-moving electrons 
in the wires. As an application of the general solution of these equations, we 
consider different kinds of electrical contacts to the double-wire system and 
calculate the direct resistance, the transresistance, in the presence of tunnel- 
ing and Coulomb drag, and the tunneling resistance. If L is smaller than the 
backscattering length Ip, both the tunneling and the drag lead to a negative 
transresistance, while in the diffusive regime {L ^ Ip) the tunneling opposes 
the drag and leads to a positive transresistance. If L is smaller than the 
phase-breaking length, the tunneling leads to interference oscillations of the 
resistances that are damped exponentially with L. 
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I. INTRODUCTION 



One-dimensional (ID) electron systems, such as they occur in semiconductor quan- 
tum wires, are in the forefront of research in modern condensed-matter physics. In 
submicrometer-long quantum wires at low temperatures the electron transport occurs in the 
ballistic regime^ and the wire conductance reaches its fundamental value of Go = c^/tt^. On 
the other hand, in sufficiently long wires the conductance is limited by scattering processes. 
If quantum-interference effects are neglected, as is the case when the inelastic scattering 
dominates, the conductance is given by cr/L, where L is the wire length and cr the conduc- 
tivity described by the Drude expression a — e^riTtr/m, where n is the electron density, Ttr 
the transport time, and m the effective mass of the electron. This regime is referred to as 
the diffusive transport. 

Modern technology allows to create various systems comprising two quantum wires put 
closely to each other so that the tunneling of electrons between the wires and/or interlayer 
electron-electron interaction is essential. Both these effects give rise to coupling between 
the electron sub-systems in single wires and in that way modify their electronic properties. 
This renders the coupled double- wire systems a subject of interest. In the past years, 
there have been experimental and theoretical studies of ID-ID tunneling^'^ and electron 
transport^"^^ along the wires of such systems. Investigations of the transport are mostly 
devoted to interlayer tunneling in the purely ballistic regime and in connection with the idea 
of electron- wave coupler^'^. On the other hand, there are theoretical papers^^^^^ describing 
the momentum transfer between the wires due to interlayer Coulomb interaction and the 
corresponding interlayer transresistance (Coulomb drag). Calculations of the Coulomb drag 
have been done both for the diffusive^^~^° and ballistic^^ transport regimes, as well as for 
the regime in which the electron sub-systems are described by the Luttinger liquid modeP^. 

Despite this progress, there is a substantial lack of description of the electron transport 
in coupled quantum wires. Even if we accept the concept that the electrons are described 
by a normal Fermi liquid, two important questions arise. The first one is how to describe 
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the electrical properties when both tunneling and interactions of the electrons (with each 
other and with impurities or phonons) are essential. The second one is how to bridge the 
gap between the ballistic and diffusive transport regimes in such a description. 

In this paper we present a linear-response theory of electron transport in coupled quan- 
tum wires that gives a reasonable answer to both questions stated above. We consider two 
parallel, tunnel-coupled ID systems of degenerate electron gases adiabatically contacted to 
four equilibrium reservoirs, as shown and labeled in Fig. 1. This general scheme of a four- 
terminal device may describe both planar^'^'^'^^ and vertically coupled^'^'^^'^^ double-wire 
devices. We take into account the interaction of electrons with themselves as well as that 
with acoustical phonons. We start from the quantum kinetic equation and finally transform 
it to a set of linear differential equations describing the distributions of the local chemical 
potentials for the systems of forward- and backward-moving electrons along the wires. The 
boundary conditions for such equations arc dictated by the Landauer-Biittiker-lmry theory. 
This transformation is justified from microscopic calculations, which also give us expres- 
sions for the characteristic times, associated with the interactions involved, that enter the 
equations for chemical potentials. 

As an application of our transport theeory, we analyze in detail different kinds of electrical 
contacts to the double-wire system. First we consider the case when the voltage is apphed 
between the ends of one wire and calculate the "direct" resistance of this wire as it is 
affected by the presence of the other one, as well as the transresistance, i.e., the resistance 
associated with the voltage induced at the ends of the uncontacted wire. Details about 
experimental measurements of the transresistance in such systems can be found in Rcf. 23. 
Next we consider the case corresponding to the tuneling measurements^, when the voltage 
is applied between the wires, and calculate the tunneling resistance. Several previously 
obtained theoretical results for such quantities (some of them are for coupled 2D systems) 
follow from our theory as limiting cases. A brief account of the main results appeared in 
Ref. 23. 

The paper is organized as follows. In Sec. II we present the microscopic model and derive 
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the equations for the local chemical potentials in the layers. In Sec. Ill we solve these equa- 
tions, apply the obtained results to the calculation of the direct resistance, transresistance, 
and tunneling resistance of the double- wire system, and describe possible transport regimes. 
Concluding remarks and discussion of the approximations made are given in Sec. IV. The 
Appendix contains detailed microscopic calculations and expressions of the characteristic 
times entering the equations for chemical potentials. 

II. FROM QUANTUM KINETIC THEORY TO LOCAL DESCRIPTION 

Consider two homogeneous ID quantum layers of length L, labeled left (/) and right (r) 
along the x axis, see Fig. 1. The quantum kinetic equation for the density matrix p reads 

dp/dt + {i/n) [Ho + Hc + He-ph, p] = 0. (1) 

Here we assume that electrons interact with each other via the Coulomb field He and 
with acoustical phonons He-ph- Elastic scattering is neglected, i.e., we assume ideal wires. 
The unperturbed Hamiltonian includes both the kinetic and potential energy operators. 
Below we use the basis of the isolated I and r layer states Fi{y, z) and Fr{y, z) and assume 
that only the lowest level is occupied in each layer. In this basis the potential energy is the 
matrix 

h^{A/2)a, + Ta,. (2) 

Here ai are the Pauli matrices, A is the level-splitting energy, and T the tunneling matrix 
element characterizing the strength of the tunnel coupling. Such tight-binding description 
is often used in application to two-level systems. 

The kinetic equation can be written^^ as one for the Keldysh's Green's function G~~^. 
Below we consider the case when the characteristic spatial scale of the electronic distribution 
is large in comparison to the Fermi wavelength nh/pF and use the Keldysh's matrix Green's 
function in the Wigner representation G~^^ {p, x) , where p and e are the momentum and 
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energy and t the time. The time dependence of G~~^ is not essential in the following, since 
we study a time- averaged, steady-state problem. The linear response theory uses a Green's 
function of the form 



Gr{p,x) = Gp-\p)+ scrip, x), 

where a and (3 are -|- or — . The unperturbed part G^^°'^(p) is given by 

cP-^ip) = miGiip) - Gf (p)], gP^-{p) = (/(£) - 1) [G^ip) - Gf (p)]. 



(3) 



(4) 



Gp--{p) = G^ip) + G^^ip), &r\p) = -cm + G'-r^ipl 



(0)++, 



V(0)- 



(5) 



where f{e) = [1 -|- e^^ A')/^B^e] jg the equilibrium Fermi distribution function and G^'^ are 
the retarded and advanced Green's functions which satisfy the equations 



(6) 



here ts'^ip) are the self-energy functions. 
The linearized kinetic equation reads 



^l^P^-^SGs + {p,x) \ + i 



h,6Gs^{p,x) 



h j d ^ d ^(0)-+, ' 



^,gP ^{p) =t6I{6,p,x). 



(7) 



Here {...} denotes anticommutators, Vp = PiVip + Pr-Vrp is the diagonal matrix of the group 
velocities, and Pi = {1 + d^j2 and = (1 ~ ^z)!"^ are the projection matrices. In this 
paper we consider the case of equal group velocities in the layers, with vip — Vrp — Vp — p/m. 
If a magnetic field is applied perpendicular to the wire plane though, vip and Vrp become 
different. Further, (p is the matrix of the self-consistent electrostatic potential arising due to 
perturbation of the electron density. In the mean-field (Hartree) approximation this matrix 
is diagonal. Finally, the generalized collision integral X is given by^^ 
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(8) 



where all Green's functions G"^ and self-energy functions S"^ have the same arguments £, 
p, and X. This corresponds to a quasiclassical description of the scattering. However, the 
matrix structure of G"^ and T,"'^ remains important and Eq. (7) is not reduced to a classical 
Boltzmann equation. Since we consider the interaction of electrons with each other and with 
acoustical phonons, the corresponding lowest-order contributions to the self-energy are given 
by the diagrams of Fig. 2. We neglect the exchange part of the Coulomb interaction for 
the following reasons. The first-order exchange contributions do not influence the imaginary 
part of the self-energy and are not, therefore, essential for the calculation of the coUision 
integral. The second-order exchange contributions are small as compared to the second-order 
direct Coulomb contributions, represented by the diagram of Fig. 2 (b), if the momentum 
transfer q is small in comparison with the Fermi momentum. Finally, there is no exchange 
contributions to the inter layer Coulomb interaction. 

We consider low temperatures and degenerate electrons. We also assume that the Fermi 
energy is large in comparison with both the tunneling matrix element T and the level splitting 
A, thereby neglecting the difference between the electron densities in the layers. We sum 
up Eq. (7) over the electron momentum p in the regions of positive (or forward) and 
negative (— ) (or backward ) group velocities and introduce the nonequilibrium part ge{x) 
of the energy distribution function in the manner 

l^Km-iP,-). (9) 



s / dp 
9e W 



Since 6G^^ is essentially nonzero only in narrow intervals of energy and momentum near 
the equilibrium chemical potential /i and Fermi momentum pp, we can replace \vp\ in this 
equation by the Fermi velocity vp, common to both layers. The integration in the -|- and — 
regions in Eq. (7) removes the contributions proportional to the potential matrix (p{x) and 
we obtain 

±VF—g^{x) + ^ [KgH^]] = Si±{e,x) (10) 

where the collision integral SI±{e, x) — {27rh)~^ dp\vp\ 51{e,p, x) depends on both g'^ and 
g~, since it accounts for both forward and backward scattering processes. However, when 



we integrate Eq. (10) over the energy, the diagonal part of di±{s,x) vanishes for forward- 
scattering contributions, and only the backscattering contributions remain, see below. In 
contrast, the forward-scattering contributions for the nondiagonal part of the collision inte- 
gral are not eliminated by the energy integration. 

The matrix kinetic equation (10) is equivalent to eight scalar equations for the four com- 
ponents of g'^ and the four ones of g~, corresponding to forward- and backward-propagating 
electrons, respectively. These equations must be accompanied by boundary conditions con- 
necting the components of g"^ with the quasi-equilibrium distribution functions of the four 
leads which the quantum wires are contacted to, cf. Fig. 1. The distribution functions of 
the leads are defined by the four chemical potentials nu, nir, ii2h and ii2r- If we assume that 
the potentials in the contact regions are sufficiently smooth in comparison with the Fermi 
wavelength but abrupt enough as compared to the characteristic scale of the electronic dis- 
tribution, we can apply Eq. (10) in the contact region as well. It gives us the conditions of 
continuity for all components of ^£ (x) across the contact regions and we obtain 

^+(0) = -^^[PiSyiu + Pr5li,rl 9e{L) = -^^[PlSfi^l + PrSfi2r], (11) 

with S/iii — /III — II, etc. The forward- and backward-propagating states are "connected", 
respectively, to the leads 1 and 2. The nondiagonal components vanish at the contacts 
because the tunneling is absent outside the region x = [0,L]. 

The problem described by the matrix equation (10) and the boundary conditions (11) 
can be considerably simplified and solved analytically if we assume that both backscattering 
and the intcrlaycr tunneling occur much less frequently than the scattering of electrons inside 
the layers and inside the -|- or — regions. The tunneling can be made weak if, for example, 
the potential barrier between the wires is thick enough. As concerns the backscattering, this 
condition is often fulfilled at low temperatures for both the electron-electron and electron- 
phonon scattering mechanisms. In the first case, the backscattering probability contains a 
factor [Ko{2ppa)]'^ , where Kq is the modified Bessel function and a is the wire width. This 
factor is exponentially small for 2pFa > 1. The acoustic phonon-assisted backscattering 
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gives a small contribution in comparison with the electron-electron forward-scattering due 
to the smallness of the electron-phonon coupling constant. In addition, this backscattering 
is suppressed at very low temperatures Tg < 2pFS, where s is the sound velocity. If the 
stated conditions are fulfilled, the diagonal part of the energy distribution function of the 
electrons have a Fermi-like energy dependence because any quasi-equilibrium Fermi function 
locally satisfies a kinetic equation with the electron-electron and electron-phonon collision 
integrals. It means that the diagonal part of gs^{x) is given by the following equation 

[9ii^)]n = '^^f^t('^) (12) 

where 5fj,f{x) = ^f{x) — fi {j = l,r) do not depend on the energy. The quantities fif{x) 
have the direct meaning of local chemical potentials for the layers / and r. It is convenient 
to introduce also the nondiagonal components of the chemical potentials fJ^^ix) — fJ'^ix) — 
ifiy{x) and fifi{x) = f^ui^) + ^(^vi^)^ by writing the whole chemical potential matrix as 

Sfi^ix) ^ jde gt{x). (13) 

Below we drop the symbol "5" in 5fi^{x) and the contact potentials S^n^r and 5^2i,r since 
all chemical potentials arc counted from the same equilibrium value yU. 

Substituting Eq. (12) into Eq. (10) and integrating the latter over the energy, we finally 
obtain eight coupled, first-order differential equations for the eight components of fJ'fix): 

± df,^/dx + (/.f - //f )(1/Zp + 1/Id) - (/^,^ - f^^)/lD - 2tFf^^ = 0, (14) 

±dfif/dx + {tif - iJ^){l/lp + I/Id) - {fit - fif)/lD + ^tFfi^ = 0, (15) 

±df^^/dx + 5p/x± + = 0, (16) 

±di^^/dx - 6fi^^ + ii^/lc + tpiiJ^f - iJ^t) = 0. (17) 

Here tp = T/hvp and 6f = A/hvp- The boundary conditions for all potentials follow 
from Eqs. (11)-(13) and are /U^(0) = /jLu, lijiL) = ^21, A*^(0) = /^ir, A^r (-^) = A*2r and 
//+^(0) = //~„(L) = 0. The characteristic lengths Ip, Id, and Ic result from the coUision 
integral 5i±{e, x), evaluated to the lowest order with respect to the tunneling matrix element 



T, see Appendix for details. They are expressed, respectively, through the phonon-assisted 
ID transport time^^ tp, the ID Coulomb drag time^^'^° td, and the phase-breaking time tc 
describing the suppression of tunnel coherence, as Zp = 2vfTp, Id — '^vfTd, and Ic — vpTc- 
The transport time Tp is common to both layers since we assume that the confining potentials 
for the wires I and r are almost identical. The analytical expressions for the Tp, td and 
Tc are given in the Appendix. All characteristic lengths are sensitive to the temperature 
T and the level splitting A. It is essential that Ic, which is controlled by electron-electron 
interaction, is always much smaller than Id and Ip. On the other hand, depending on the 
temperature and level splitting one can have different relations between Id and Ip: both 
cases, Ip <^Id and Ip ^ Id, are possible. 

Equations (14)- (17) with the stated boundary conditions give us a complete description 
of the electrical properties of double quantum-wire systems in a wide range of regimes 
starting from the purely ballistic transport regime L <^lc to the diffusive transport regime 
L ^ Ip,Id- The local currents flowing in the layers j — l,r are expressed by 

Jj{x)^Go[fit{x)-fij{x)]/e, (18) 

and the local tunnel currents are proportional to Tii^{x). 

Below we present the general solution of Eqs. (14)- (17) and describe two important cases, 
that of long systems, with L ^ Ic, and that of short systems with L ~ Ic- To characterize 
the effects of drag and tunneling, we then consider different kinds of electrical contact to 
the double-wire system. First we consider a typical setup for the drag measurements, when 
the current J,. = Jr{0) = Jr{L) is injected in wire r (''drive wire") while no current is 
allowed to flow into wire I, Ji{0) — Ji{L) — 0, and calculate the transresistance Rtr deflned 
as Rtr — [^J^ll — lJ'2i]/eJr as well as the "direct" resistance R — [/iir — ^i2r\/&Jr- Next, 
we turn to the tunneling measurements^, when the voltage is applied between the wires. 
We consider both the symmetric setup, when all four ends of the wires are connected to 
external sources, with — fj,2i and /iir — ii2r, and the non-symmetric one, when the 
voltage is applied between the ends 11 and 2r while the remaining ends are not contacted, 
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Jr(0) = 0, Ji{L) = 0. For each of these cases we calculate the tunneling resistances Rts 
(symmetric) and Rm (non-symmetric). Both of them can be defined as [jiu — ^JL2r\/^JT-i 
where Jt is the total current injected. Jt is equal to 27^(0) and Ji(0) for the symmetric and 
non-symmetric contacts, respectively. 

III. RESULTS 

Since Eqs. (14)-(17) are hnear, their general solution is easily obtained as 

(x) = (1 + L/lp)-^ [{im + /ilr) (1 + (i^ - X)/Ip) + {H2l + l^2r)xllp\ 

± (A^e^^'^ + B+e-^^"") (19) 

i 

2//^- (x) = (1 + L/lp)-^ [{l^u + fiir){L - X)/Ip + {H2l + l^2r) (1 + x/lp)] 

± ^ (A"e^^^ + Bre-^^"") (20) 

i 

/^^(^)=E(^^^'' + A^e-^i (21) 

i 

,ii^)--^(-^.Cte^'.^.Oi.-^'). (22) 
Here Aj = -s/yl and y — fji are the solutions of the cubic equation 

-Atl[A{tp/lcf + {2/lch){lc' + Si)] = 0, (23) 

where ^ = /p^ -|- 2/^^. The coefficients ylf , 5^^, C^, and are to be found from Eqs. 
(14)-(17) and the relevant boundary conditions. Below we use the property Ic <^ Ip, Id and 
the condition of weak tunnel coupling tp <^ Iq^ to simplify this procedure. Then the three 
roots of Eq. (23) are easily obtained as 

Ai = A~2(l//2 +i//^/,)V2^ A2,3 = A±~l//c±^(^F, (24) 

where we introduced the tunneling length It = vpTr- The tunneling time tt, which contains 
a resonance dependence on the level splitting, is defined by 
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''^'^^^\' + {h/rcr ^^^^ 
The root A describes long-scale variations of the chemical potentials while X± corresponds 
to short-scale variations. Accordingly, we consider the regimes that follow. 

A. Long wires, L 2> Ic- 

This length range comprises the region from the "pseudo-ballistic" {Ic <^ L ^ Ip,Id) 
to the diffusive (L ^ Ip) regimes. All solutions containing A± exist only in short regions in 
the vicinity of the contacts. They are evanescent inside the wire region and not essential in 
the calculation of the currents. Considering only the solutions involving A, we find 

Ml^ii - liir)P{L - x)/P{L) ± {ii2i - l^2r) smh \x/P{L) (26) 

2hJXx) = (1 + L/lp)-^ [{nn + nr,){L - x)/lp + {^21 + /^2r) (1 + x/lp)] 

Hl^ii - l^lr) sinh A(L - x)/P{L) ± {ii2i - ii2r)P{.x)/P{L), (27) 

where P{x) = (1 + 2/i//t) sinh \x + A/i cosh \x. The same expressions can be obtained from 
the four coupled balance equations compactly presented as 

± d^xf/dx + {^xf - /xj)(i//p + i/Id) 

-{pip - pif,)/ln + {pif - 14.) /It = 0, (28) 

where j — l,r and / 7^ j. These equations follow from Eqs. (14)-(17) in the limit L » Ic, 

when one can neglect the derivatives d/i^^^/dx in comparison to /^^^^//c- With Eqs. (26) 
and (27) we obtain 



2 + L/lp + (1 + h/h)^^^ tanh(AL/2)] , (29) 



L/lp - (1 + h/h)^^^ tanh(AL/2)] , (30) 
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For L -C IpiIt we have R ~ uTi/e^ = Gq^, while the transresistance is given by 
Rtr ~ —{7rh/e'^)L[l/lr) + 1/21t\. As seen, Rtr is small, always negative, and propor- 
tional to the wire length L multiplied by a sum of drag and tunneling rates. If one neglects 
tunneling, the resulting expression for Rtr, with td given by Eq. (All), describes the 
Coulomb drag in the ballistic regime previously investigated^^ by Gurevich et ai. When L 
increases and the electron transport becomes diffusive (L » Zp), we obtain, for XL/2 <S 1, 
R ~ {7rh/e^)L[l/lp + 1/Id]- This resistance, if one omits the drag contribution, is expressed 
in terms of the usual Drude conductivity a — L/ R — e^lp/irh — e^riTp/m. The correspond- 
ing transresistance is 

Rtr = -^{L/Id) [l - (V^o)'l , Lo = {GlllT/lof' ■ (31) 

Expressing l^ and It through the drag transresistivity Trh/e^lo and the tunneling conduc- 
tance Gt = g^Pid/tt = le^ /-kTiIti where piD is the ID density of states at the Fermi level, 
one can see that Eq. (31) formally coinsides with that obtained in Ref. 27, where a com- 
petition of drag and tunneling effects was investigated for double quantum-well systems. 
For XL ~ 1, the transresistance is large and comparable to the direct resistance, because 
a considerable fraction of the current penetrates the I layer due to tunneling. This regime 
for double quantum wells has been investigated both experimentally^^ and theoretically^^. 
If one neglects the drag and assumes the diffusive regime (L ^ Ip) with weak tunneling 
{Ip <^ It), Eqs. (29) and (30) describe the results obtained in Ref. 29. For XL 1 we have 
Rtr — R — {'Kh/eHp){L/2) — L/2a. This is the case when the current, though injected 
only in one layer, is equally distributed among the layers due to tunneling. 

Figure 3 shows the length dependence of the transresistance calculated for different rel- 
ative contributions of the Coulomb drag and tunneling. The transresistance is negative for 
small L but always changes its sign and becomes positive as L increases and the backscat- 
tering occurs often [see also Eq. (31)]. This behavior can be explained with the help of 
the balance equation (28), which shows that the tunneling tends to decrease the difference 
between lif and while the backscattering tends to decrease the difference between yil^^ 
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and Thus, for /i;"^(0) = /X;^(0) and nt^L) = iiJ^L) the change of (fJ^f ), with x, 

is opposite to that of /i^ at small L and becomes the same as that of /i^ as L 
increases, leading to the change of Rtr — [A*f (0) — ixf{L)]/eJr from negative to positive. 
This transition occurs at smaller L/lp if the tunneling is stronger (larger Ip/lr) and the 
drag weaker (smaller Ip/lo)- 

Although It is normally longer than Zp, the opposite condition can also be realized. 
A particularly interesting transport regime, corresponding to long quantum wires with- 
out backscattering, occurs in tunnel-coupled magnetic edge states^°'^^, since an edge state 
represents a ID system where the electrons can move only in one direction. Assuming 
l/lp = I /Id = in Eqs. (29) and (30), we obtain the result of Ref. 31 in the form 

TTn TTu 

= — [l + (l/2)tanh(L//r)], i^rii = tanh(L//r). (32) 

Consider now the behavior of the tunneling resistances. With Eqs. (26) and (27) we 
obtain 



1 + (1 + h/hf^^ coth(AL/2)] (33) 



and 



2 + L//p+ (l + /T//i)^^^coth(AL/2) , (34) 



for symmetric and non-symmetric contacts, respectively. For conditions XL/2 ^ 1 we have 
Rts — Rtu — (7^^/2e^)(/T/-f'), i-e., the tunneling resistances depend only on the ratio of the 
tunneling length to the wire length. This is because the regime of AL/2 <S 1 corresponds 
to weak tunneling and the chemical potentials f^^{x) and lif{x) arc close to fiu and ii2r, 
respectively. With the use of the tunneling conductance Gt (see above) one can rewrite the 
expression for the tunneling resistances in a more transparent way: Rts — Rtu — {GtL)'^. 
For AL/2 ~ 1, when the coordinate dependence of the chemical potentials in the layers 
is essential, Rtu is different from Rts and both of them depend on the scattering length 
Ip. The drag effect is not so important as for the transresistance: the tunneling resistances 
depend on only iilo is comparable to or smaller than both Ip and It- 
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Figure 4 shows the length dependence of the tunnehng resistances Rts and Rtu, as 
given by Eqs. (33) and (34), for several different values of the ratio Ip/It describing the 
strength of the tunneling with respect to the backscattering. The drag effect is neglected, 
1/Id = 0. As the wire length becomes larger than the backscattering length, the 1/L 
decrease of the tunneling resistance changes to either a L-independent behavior (for Rts) 
or to a linear increase (for Rtu)- In the first case the dependence on L disappears because 
all tunneling occurs near the ends. In contrast, for non-symmetric contacts the resistance 
Rtu is determined by the Ohmic resistance of the wires instead of the tunneling effects, 
and increases linearly with L. A similar effect, in applications to coupled quantum wells, is 
discussed in Ref. 29. 

B. Short wires, L <^ Ip, Id 

This length range comprises the region from the purely baUistic {I <^ Ic) to the "pseudo- 
ballistic" {Ic <^ L <^ Ip, Id) regimes. Since the electrons pass along the wires almost without 
backscattering, R is close to 7r^/e^, and Rtr is small. However, for L ~ an electron 
tunneling between the layers does not lose its phase memory completely and tunnel coherence 
effects can manifest themselves on such short lengths giving additional contributions to the 
transresistance Rtr and the tunneling resistances Rts and Rtu, accordingly the expressions 
for these quantities obtained in the previous subsection ior L <^ Ip,Id should be modified. 

A convenient analytical approach to the problem in this regime is to solve Eqs. (14)- (17) 
by iterations taking i^i'^ri^) — lJ>ii,r, fJ'^ri^) — l^2i,r and lJ>u,v{^) = as an initial approximation. 
Another way is to use Eqs. (19)-(22) directly. We obtain 

2^,+ (x) = {/In + jJiir) (1 - x/lp) {lJi2l + A*2r) x/lp 

± il^n - l^ir) [1 - x/h - 2x/It + 2{Ic/It)^{x)] ± (/^^ - l^2r) x/h (35) 

2/X^"^(L -x) = {lJ,2l + A*2r) (1 - x/lp) + {nil + nir) x/lp 
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± il^2i - l^2r) [1 - x/h - 2x/It + 2{Ic/ItMx)] ± iiiii - iiir) x/h (36) 

where 

= [l/{lc'' + 51)] [2(5p/Zc)e-"/'^ si^iSpx) + {Ic" - 4) (l - e'"/'^ cos((5i.x))] , (37) 
is an oscillating function of the coordinate x and = A/hvf. Now and i^^ij are given, 



respectively, by 



flTT 

R^—[l + L/lp + L/Id + L/21t - {Ic/21tML)\ , (38) 



and 



??7r 

i^TR = — [-L/Id - L/21t + {lc/2lTmL)] , (39) 

The contribution to i?TR coming from the term proportional to $(!/) is not small for L ~ Ic- 
It describes oscillations damped due to the factor exp(— L/Zc)- The periodic behavior can 
be described as a result of the interference of electron waves of the left and right layers along 
the length L: due to a finite level splitting A these waves have different phase velocities. 
Similar interference effects occur in the tunneling resistances 

Rrs ^ RTa ^ ^ih/L) [1 - {lc/m{L)]-' . (40) 

Prom Eq. (37) for $(-L) we see that both the transresistance and the tunneling resistances, 
being functions of AL/hvp, oscillate with level separation A. The oscillations are damped 
when the wire length L exceeds Ic so that the tunnel coherence over the wire length is 
suppressed. 

Changing A by applying a voltage across the wires would lead to oscillations of R, Rtr, 
Rts, and Rtu- Another way to change A is to apply a magnetic field B perpendicular to the 
plane of the wires^^'^^. For sufficiently weak B the results presented so far still hold with the 
phase 6fL having an additional contribution 27r0/0o, where (f)o = h/e is the magnetic flux 
quantum and (p — BwL the flux enclosed by the area between the wires. Though the double- 
wire system does not form a closed current loop, this should lead to Aharonov-Bohm-type 
oscillations in the resistances deflned by Eqs. (38)-(40). 
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In very short wires, with L -C Ici 5p^i Eqs. (39) and (40) become 

RTR^-^-^{L/ln + Lhl/2), Rts^Rtu^^L-H-/. (41) 

In this regime only a small fraction of the electronic wave packet is coherently transmitted 
from one wire to another. The tunneling contribution to Rtr follows a dependence, 
instead of the linear dependence occuring for Ic <C L <^ Ipjo, when the tunneling is 
non-coherent. The length dependence of the tunnehng resistance follows a law. 

In the investigation of the purely baUistic regime (L -C Ic) we can neglect the coUision 
integral in Eq. (10) and need not make the assumption about the smallness of the tunneling 
matrix element which was essential for evaluation of the scattering-induced contributions in 
Eqs. (14)- (17). The electron transport in coupled quantum wires in this regime is pertinent 
to the problem of electron- wave directional couplers. Theoretical studies of this problem^^-*^^, 
although rather extensive, included only a quantum-mechanical calculation of the electronic 
transmission. Below we show how the essential results of these studies can be obtained in a 
simple way from the quantum- kinetic analysis. Integrating Eq. (10), with SI±{e, x) = 0, over 
the energy and taking Eq. (13) into account, we find that the distribution of the chemical 
potentials is again described by Eqs. (14)-(17) without the terms containing the scattering 
lengths Ic, Ip, and Id- Since there is no backscattering, the solutions for /^^^^ and /i^j. are 
decoupled 

(x) = nii^r =F ifJ'ii - Aiir)r sm'^{ATx/2hvF), (42) 
I^TA^) = l^2i,r T {l^2i - H2r)rsm'^{AT{L - x)/2hvF)] (43) 

here Ay = (A^ + 4^2)^/2 j. ^ /^T'^/A'^. Equations (42) and (43) describe oscillations 
of the electronic wave packets between the layers due to coherent tunneling. A complete 
transfer of the wave packet can be achieved for A = 0. One can calculate the resistance and 
transresistance as 

R=^\l- (r/2)sin2^]/[l-rsin2^], (44) 
Rtr^-tt^t sit? V'/ [1 - r sin^ V] (45) 
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and the tunneling resistances as 

^ ^ Tr/i Trh 19, / , „, 

Rts ^^Tn-^^ ^r-^ sin-2 ijj, (46) 

Here ip = ATL/2hvF- The oscillations of these quantities occur in a way similar to the one 
described by Eqs. (38)-(40): SpL coinsides with 2ip if one replaces A by At- However, 
since the tunnel couphng is strong, the oscillations described by Eqs. (44)- (46) have large 
amplitudes. In particular, when A is small (r ~ 1), all the quantities given by Eqs. (44)- (46) 
show giant oscillations with amplitude large in comparison to Gq^. 



IV. CONCLUSIONS 

In this paper we carried out a theoretical study of electron transport in parallel ID 
layers coupled by tunneling and Coulomb interaction and contacted, at their ends, to quasi- 
equilibrium reservoirs. A linear-response, steady-state regime has been investigated, and 
the wires were assumed to be ideal, i.e., without defects and, therefore, the elastic scattering 
of electrons by them was neglected. As the most important result of our study, we found 
that a full quantum-kinetic description of the problem is reduced, with physically reasonable 
assumptions, to a set of linear, first-order differential equations describing the distribution 
of local chemical potentials for forward- and backward-moving electrons. The boundary 
conditions for the chemical potentials are determined by the potentials of the reservoirs 
controlled by applied voltages. The solution of this set was obtained analytically and allowed 
us to describe the local currents flowing in each layer from the pure ballistic regime, when 
the electrons do not suffer any scattering along the wires, to the diffusive regime, when the 
electrons experience many backscattering events during the transport. 

In particular, we applied our approach to the description of the resistance R, transre- 
sistance Rtr, and tunneling resistances Rts and Rtu of double quantum wires. The most 
important result is that Rtr, which is caused by both tunneling and Coulomb drag effects, 
depends on the wire length L non-monotonically and always changes its sign as L increases. 



17 



because in shorter wires, when backscattering is rare, the tunnehng, as well as the drag, 
leads to a negative Rtr, while in longer wires, when the transport becomes diffusive, the 
tunneling leads to a positive Rtr, in a way similar to that for coupled two-dimensional (2D) 
systems^®'^^ , and overcomes the drag as the length L increases. This sign inversion is qual- 
itatively understood from an analysis of the balance equation (28) and is mathematically 
described by Eq. (30). In the diffusive limit and for It ^ Ip, Eqs. (29) and (30) formally 
coinside with those obtained previously^^'^^ for coupled 2D systems. Besides, some recent 
studies of transport phenomena in coupled ID layers, namely transport without backscat- 
tering in tunnel-coupled edge states^°'^^ and Coulomb drag between quantum wires in the 
ballistic regime^^, constitute limiting cases of the more general results given by Eqs. (29) 
and (30). 

One should stress the importance of the phase-breaking processes that suppress the 
tunnel coherence. In our model, i.e., without elastic scattering, these processes arc proved to 
be much more frequent than the backscattering processes. This allowed us to distinguish two 
transport regimes: the pure ballistic regime, without any scattering, and the pseudo-ballistic 
one, without backscattering but with essential forward-scattering due to electron-electron 
interaction, and with suppressed coherence. If one considers just a single wire, there is no 
difference between these regimes as concerns the electrical properties: the resistance is equal 
to Gq^ in both cases. However, the electrical properties of a tunnel-coupled double-wire 
system behave differently as one passes from one regime to another, because the contribution 
of the tunnehng to the electrical properties becomes different. In the ballistic regime, as 
well as in the transition region between the two regimes (L ~ /c), all calculated resistances 
oscillate with L and with the level splitting energy A due to interference of the electron waves. 
One can vary the level splitting energy by applying either a transverse voltage across the 
wires or a magnetic field perpendicular to the wire plane. In the latter case the oscillations 
show a Aharonov-Bohm periodicity associated with the magnetic flux penetrating through 
the area Lw between the wires. The oscillations become exponentially damped as the ratio 
L/lc increases. If the tunnel coupling is strong, the oscillations have large amplitudes, 

18 



which, from a theoretical point of view, can be much larger than the resistance quantum 
Gq^. So far the experimentally observed^'^^ resistance oscillations in tunnel-coupled baUistic 
quantum wires were of small (~ 0.5 KQ) amplitude. This is not surprising because there are 
many factors which compete against the tunnel coherence. Apart from inelastic scattering 
considered in this paper, there are elastic scattering and long-scale inhomogeneities of the 
wires which would lead to a coordinate dependence of the level splitting A. If these variations 
of A are larger than the tunneling matrix element, the coherence would be considerably 
suppressed. 

We now discuss the approximations made in this paper. The main approximation is the 
neglect of elastic scattering. Since this scattering tends to be dominant at low temperatures, 
the presence of impurities in the ID channels will considerably modify the transport. The 
elastic scattering will lead to an increase of backscattering and interference between the 
forward- and backward- moving electron waves. As a consequence, R, Rtr, Rts and Rtu will 
depend on the spatial positions of the impurities in the channels and the regular expressions 
obtained in this paper will not be valid. A further development of the transport theory for 
tunnel-coupled wires in the presence of elastic scattering is therefore desirable. On the other 
hand, advances in the technology of nanostructures, in particular selective doping, can make 
it possible to achieve structures where the elastic scattering in ID channels is minimized for 
wire lengths smaller than a few microns which is the current standard of the impurity mean 
free path at low temperatures. 

Another approximation concerns the transition from the quantum kinetic equation (1) 
to the semi-classical description given by Eq. (7). It is valid when the spatial scale of the 
electronic distribution is large in comparison to the Fermi wavelength irh/pp- We have seen 
that this scale is determined either by A, given by Eq. (24) for long wires, or by A± for short 
wires. In the case of strong tunnel coupling the characteristic scale is given by hvp/Aj'- 
Therefore, the necessary requirement is fulfilled if the tunneling matrix element T, level 
splitting A, and the energy H/tc, associated with the smallest scattering time tc, are small 
in comparison to the Fermi energy. These conditions have been assumed througout the 
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paper. This also allowed us to neglect the difference between the electron densities in the 
layers and characterize the electrons in different layers by the same Fermi velocity \vp\ ~ vp- 

The assumption about the adiabatic connection of the wires to the leads, which allowed 
us to neglect elastic scattering of electrons near the ends of the wires, implies that the Fermi 
wavelength ttH/pf must be small in comparison to the contact lengths, i.e., to the lengths of 
transition from the leads to the wires. On the other hand, the oscillations associated with 
the tunnel coherence, cf. Sec. Ill B, can be seen if the contact lengths are smaller than both 
Ic and hvp/^T- In principle, both requirements can be fulfilled. 

The next approximation, which allowed us to solve the kinetic equation analytically in 
the whole range of regimes from ballistic to diffusive, is equivalent to the following state- 
ment. In each layer the forward- and backward-moving electrons can be described as weakly 
coupled sub-systems characterized by their own local chemical potentials. This statement 
is obvious for the case of pure ballistic or pseudo-ballistic transport, when these poten- 
tials are merely dictated by the reservoirs (leads) and do not change with coordinate x. 
When the backscattering becomes essential, this statement is still true if we assume that 
the forward-scattering events are much more frequent than the backscattering and tunneling 
events. For example, it is always true for magnetic edge states, where one can completely 
neglect backscattering, and the introduction of local chemical potentials (see Ref. 31) is 
well-justified. In our case, a consideration of the electron-electron collision integral allowed 
us to estimate the characteristic time of the Coulomb-assisted forward scattering, and we 
find that it is of the order of tc, which is small in comparison with both backscattering 
times Tp and td- Thus, the electron-electron interaction provides an effective mechanism 
for forward scattering and can maintain quasi-equilibrium Fermi distribution functions for 
forward- and backward-moving electron sub-systems. However, these conditions may be 
violated when the conducting channels contain impurities with short-range potentials and 
the elastic backscattering becomes important. 

Our evaluation of the characteristic scattering times from the collision integral has em- 
ployed only the lowest-order essential contributions of the electron-phonon and electron- 
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electron interactions, given by the diagrams of Fig. 2 and leading to collision integrals with 
scattering amplitudes in the Born approximation. While it is normally^^ good for electron- 
phonon interaction due to the weakness of the coupling constant, a rigorous evaluation of 
the electron-electron part requires also a consideration of higher-order contributions, given 
by more complex diagrams, because the ratio of the Bohr energy to Fermi energy Sp, which 
is the parameter of the perturbation expansion for the Coulomb interaction, is not small. 
Nevertheless, using the Born approximation in the evaluation of the drag time is still rea- 
sonable if the momentum 2pF transferred in backscattering is large and the electron-electron 
backscattering probability is small. As concerns r^, it is determined by forward-scattering 
processes with small momentum transfer and the Born approximation is not well justified^^. 
On the other hand, our theory leads to a non-divergent expression (A14) for tc and gives, 
for typical parameters of the electron system, physically reasonable values. We remind that 
in our theory both ksT^ and A are much smaller than ep- Therefore, one may expect that 
Eq. (A14) provides a correct order-of-magnitude estimate of the phase-breaking time caused 
by electron-electron interaction. 

Finally, we stress that the results obtained in this paper hold for a normal Fermi-liquid 
state of the electron system. If the electrons in the wires are in the Luttinger-liquid state^*^, 
these results have to be reconsidered. 
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APPENDIX 

Below we give a microscopic calculation of the characteristic times Tp, td, and tc- The 
coordinate index x in the Green's functions and self-energies is omitted and Ti is set equal 
to 1. The normahzation lengths are also set equal to 1. The electron-phonon self-energies 
given by the diagrams shown in Fig. 2 (a) are explicitly expressed as 
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^ikeiP) = ^ E / Z^Tf,e -JP- ^D!f(^^ Q)M;r\Q)M;r/\-Q), (Al) 
Q 

where the unperturbed Green's functions of phonons D~~^ and D~^~ (we do not need D 
and in the following) are given as 

iDp{uj, Q) = 27r [Nq5{uj ^ sQ) + (1 + Nq)5{uj ± sQ)] {A2) 
and the matrix elements of the electron-phonon interaction are 

Jr'i%^^^) =J J dydzFf{y,z)e'''yy+"^^\ {A3) 

We use the expression sQ for the phonon energy, where Q — |Q|, Q — [qAyilz] is the 
phonon wave vector and s the velocity of sound. Further, Nq = 1 / [exp{sQ / kBTg) — 1] is the 
Planck distribution function, p is the material density, and Ei is the deformation potential 
constant. 

The electron-electron self-energies given by the diagrams shown in Fig. 2 (b) are ex- 
pressed as 

jij'i P',Q 

where / = for a = /3 and / = 1 for a ^ P, the factor of 2 comes from the spin summation 
in the "loop" . Here 

Mp(g) = (2eV6) III Jdydy'dzdz'Ko{\q\\r-r'\)Ff{y,z)F}{y',z') (A5) 

are the matrix elements for electron-electron interaction, e is the dielectric constant, Kq is 
the modified Bessel function, and |r — r'| = [{y — y'Y -\- {z — z'Y^I'^ . 
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For the evaluation of the colhsion integral we use Eqs. (3)-(6) and express the non- 
equilibrium part of the matrix Green's functions according to [see also Eq. (9)] 



SGf{±\p\) = giGiip) - G§{p)gi. {A6) 

The collision integrals are evaluated below assuming of weak tunneling, when the non- 
diagonal contributions of G£'^{p) are neglected. This approximation is valid when the 
tunneling matrix element T is small in comparison to the imaginary part of the self-energies 
and when the level splitting |A| is small in comparison to the Fermi energy. Both require- 
ments are assumed fulfilled. Then the components [g£']jj' enter only in the corresponding 
parts S[I{e)]jj' of the collision integral. 

Now we calculate the diagonal parts for the electron-phonon scattering contribu- 
tion to the collision integral. Taking the self-energy given by Eqs. (A1)-(A3) we find 
/ deS[i^^'^{e)]jj — —{lJ>f — iij)/2Tpj where j — l,r. The phonon-assisted transport time is 
given by 



P^Te ' " P,.(P>t;-,<o) " J 4 sinh^(.g/2Te) 



X 



{A7) 



where we defined G'^ = G^ — G^. For further evaluation of Tpj we use the free-particle 
(unperturbed) Green's functions, i.e., G^'i{p) = [s - A/2-p^/2m ± iO]-^ and G^;^{p) = 
[e + A/2- p^/2m ± iO]-\ and obtain 



-1 _ El 2 r A + '"Fqi/2 - sQ/2) - /(/x + VFqi/2 + sQ/2)] 

- psksn h ' " 4 sinh^(.Q/2fc^r.) 

(AS) 

where qi — q — 2pF is a small variable. In the limit ksTg ^ ms^ the integral over qi is 
easily calculated. In addition, if ksTe s[(7r/a)^ -|- (2pi?)^]^/^, where a is the wire width, 
the scattering becomes quasi-elastic, and Eq. (A8) is reduced to a known^^ result 
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_i _ 2ElkBT, 



J J dydzFf{y,z). (A9) 



Since we assume that the confining potentials for the layers are almost identical, the dif- 
ference between Tpi and Tpr is neglected: Tpi — Tpi — Tp. A numerical estimate, using Eq. 
(A9) and GaAs material parameters, gives Tp^ ~ 10~^ ksTf,. 

The electron-electron scattering contribution to the diagonal parts of the collision inte- 
grals gives the Coulomb drag terms: / de5[Ij^'^{e)]jj — —{jJ^f — ijJ)/2td + [jJ^fi — IjlJ,)/2td, 
where / ^ j. The drag time is given by 



^ p,p ,q 



The sum here must be evaluated forp>0,p — g'<0, p'<0, and p' + q > 0. The evaluation 
of Eq. (AlO) using the free-particle Green's functions gives a simple result 

_i knT {A/2kBT,f 
nvp smh [A/2kBTe) 

One can estimate Mf~^{2pp) as {2e'^ /e)Ko{2pFw), where w is the distance between the 
centers of the wires. If 2pFW -C 1, which is easily achieved for w ~ 30 nm, Ko{2pFw) is 
exponentially small. 

Finally, we calculate the electron-electron part of the nondiagonal components of the col- 
lision integral. Since the main contribution to it comes from the forward-scattering processes 
{\q\ <S Pf), only such processes are considered below. The integral of [I±'^{s)]jj' {j 7^ /) 
over the energy e can be reduced to a sum of three terms characterized by three different 
statistical factors: 



/ deS[it\e)],r = -2/ de[gi],,, ^ ^p J J 



de'duo 
(27r)2 



py,?(p>o,p-?>o) 

X { KT"" + Kf") - f{e' + a;)] + f{e - u:)[f{e' + u) - f{e')]] 
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+ (A#^^ + A2^^^) f{e' + f{s')] + {A^,fr + Af^^^) f{s - u)[f{s' + f{e')]} . 

{A12) 

In Eq. (A12) we used the shortcuts 



X 



{A13) 



and neglected the terms with a = (3 and 7 = 5 because they vanish after the summations over 
p and p', respectively. Since [gijir = [gsVri = 9u,s - wt,e, one can see that / de5[i^''{e)\jj> = 
!de5[il-^{e)]],^. 

Calculating the integrals in the expression (A12) within the approximation of the free- 
particle Green's functions, we find that the third term on the right-hand side of Eq. (A12) 
vanishes. The first term diverges for A = but it is close to zero for A 7^ and can 
be neglected in the following. In contrast, the second term gives a regular contribution, 
which can be represented, on account of Eq. (13), as / de5[I^^{e)\jji — —nfji/rc- The 
"nondiagonal" relaxation time tc (we take into account only its real part) is given by 

e^S'^A A 

where 



J J J dydy'dzdz'ln |r - r'\F^{y, z)[Ff(y', z') - F^{y', z')]. (Alb) 

In the calculation we took into account Mff^{q) ~ M^~^{q) and qa <^ 1. The last property 
allowed us to use the approximation Ko{x) ~ —[C + ln(a;/2)], where C is Euler's constant; 
we found Mif^^q) — Mf~^[q) ~ {2e^/e)S. The overlap integral S can be approximated, to 
a good accuracy, by ln{w/a). 

If A » A^ksTe., the relaxation rate Tq^ given by Eq. (A14) is temperature-independent 
and proportional to |A|. For A <^ AksT^^ Tq^ is proportional to Tg. A comparison of Eq. 
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(All) and Eq. (A14) shows that rc is always much smaller than td, since tc is controlled by 
forward-scattering processes and does not contain the smallness associated with the factor 
[^0(2^^^)]^. A numerical estimate also shows that tc <S tp, because of the weakness of 
the electron coupling to acoustical phonons. For this reason we neglected the contribution 
of electron-phonon scattering to the nondiagonal part of the collision integral. 
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FIGURES 



FIG. 1. Schematic representation of two coupled quantum wires. 

FIG. 2. Feynman diagrams describing the contributions of electron-phonon (a) and elec- 
tron-electron (b) interaction to the self-energies. 

FIG. 3. Transresistance Rtr as a function of the wire length L. The solid curves correspond 

to Ip/lo = 0.1 (weaker drag) and the dashed ones to Ip/lo = 1 (stronger drag). Each curve is 
marked by the value of Ip/lr- 

FIG. 4. Tunneling resistances for symmetric (solid) and non-symmetric (dashed) setups as a 
function of the wire length L, at l/ln = 0. Each curve is marked by the value of Ip/lx- The inset 
shows the currents (arrows) injected in and coming out of the wires (broad lines) for both cases. 
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Fig. 1 . O. E. Raichev and P. Vasilopoulos, Electron transport 




Fig. 2. O. E. Raichev and P. Vasilopoulos, Electron transport.. 



